481 

Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 
nCIDES-Iin. Editor: G.S. Dulikravich. Washington D.C.. October 23-22. 199L 


r 


OPTIMAL INTERACTION OF INDENTER WITH INHOMOGENEOUS PLATE 

Valery N. Aptukov N92-13964 


1 


Institute of Continuous Media Mechanics 
UB of the Academy of Sciences USSR, 
Koroleva st. , 1, Perm, 614061, USSR 

ABSTRACT /> ' 



Consideration is given to a new class of problems dealing 
with an optimal design of inhomogeneous plate during dynamic pene- 
tration of the rigid indenter. The quality criterion of the pro- 
cess is defined by the specific mass of tne target, which absorbs 
the given kinetic energy of the indenter. Parameters of control 
are expressed in terms of mechanical characteristics, i.e. distri- 
bution of density p and the related hardness H across the plate 
thickness. The maximum principle of Pontryagin are used to search 
for piece-wise continuous control function. With consideration of 
impact conditions and characteristics for a given class of mate- 
rial an optimal target structure criterion has been estimated for 
engineering application. 


INTRODUCTION 


The problem of searching for mechanical characteristics of 
inhomogeneous plate subject to impact of a rigid body has been 
stated first in til in The framework of theory of optimal control. 
This study employs Pontryagin principle of maximum 1 23 to obtain 
an optimal structure for a plate with minimal thickness and 
prescribed specific mass. At present a considerable attention is 
focused on the problem of structure optimization as applied to the 
case of inhomogeneous plate of a minimal specific mass using both 
linear HCp 0 - Ap + B [3] and nonlinear HCpj - ipCpO [4] relations. 
An approximate approach to the analysis of penetration process, 
based on the empirical relation [53 allows to obtain rather simple 
criteria for structure optimization. 

ANALYSIS 


1. Formulation of the problem. 

According to the applied theory of the plate specific resis- 
tance p , penetration of the rigid indenter can be expressed [53 

clS 

p = H + kpu 2 C 1 . 15 

where H is dynamic hardness; k is the shape factor of the indenter 
head (in case of a tapered head k= sin 2 a, a is a half angle of the 
cone opening); v is tne current penetration rate. 

Tne equation of motion for indenter is given as [33 

U/20M dCv z 0/dL = -2n J pCx,L0 rCPdrCfO/d{ d£, Cl. 2) 

O 

•where L is the current penetration depth; r(£) is the expression . 
for generating line of axisymmetrical indenter; £ = L-x is the J 
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Coordinate relative to the tip of the indenter Cfig.l). 

Distributions of density and hardness are assumed to meet 
following requirements 

pQQ, Q = < pCxi>: p t < pCx 0 < p g , CxQ [0,L ]} > 


the 


WS*, 3? = { HCxd: H < HCx 0 < H 


Cx c [0,L k D ) Cl. 3) 


We shall restrict ourselves to a class of materials for which 
there exists a one-to-one mapping <p of the set Q in the set 3E 
H - ipCpDy p t H i , p z $ H z with dip/dp > 0. According to the quality 


criterion stated below, materials inconsistent with the latter 
conditions are considered inadequate. Further it is assumed that 
each of the materials is plastic enough and impact velocities lie 
within the range, in which application of the relation Cl. 2) 
proved to be valid. 

The boundary conditions for Cl. 2) imply that the indenter 
moves with the initial penetration rate vl h=0) = v and reaches 

some unknown finite penetration depth L k for which uCL k ) = 0. 

We shall concentrate on the case with the plate thickness being 
equal to the finite penetration depth b = L k . 

Specific mass of such target is taken as the principle crite- 
rion of quality L k 

J = min ( f pCx ) dx } Cl. 4) 

pcQ l o / 


1 


2. Conical indenter. Linear relation H = <pCp 0. 

Let us consider the technique of applying the maximum prin- 
ciple [23 to a number of particular problems. Within the framework 
of theory of optimal control the problem may be expressed in the 
form of C_ <y‘ = v , y, y 3 ) is the vector of phase coordinates, t 

U/OB t z + CA + My 1 ) y 2 ], 

= y 3 , dy 3 /dt = p, C2.1) 

reasonable to introduce additional 
phase coordinates y, y , since the right side of the equation of 
motion involves an explicit form of functional Cl. 3), leading to 
condition dyj Q /dt- -dh/dy 0 / 0. 

The mapping H = pCp) is assumed in the form of linear appro- 
ximation 

HCpd = ApCxd + B t 

A = CH z - N^/Cp m - Pi )y 6* CH i p z - H zPi d/Cp z - p^y C 2.2) 

At the initial moment of time the vector of phase coordinates 
remains fixed t = 0 : y 1 = y^ , y 2 = y 3 = 0. The finite vector value 

space 


C2.3) 

van- 


two-dimensional variety S, Euclidean 


y k belongs to a smooth, 
with dimension n = 3 

S k : F k Cy k' y k' E K = 0 

.The condition of transversal ity for the vector of conjugate 
[ables yj yields two relations 


= L 


is une Lime analog j 
dy'/dt = -E [ 

dy z /dt 

where E = C4n/M) Ctgcu 2 . 

In the followinq it is 
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t = t k : y z = 0, w 3 = 0 

The equation for the conjugate variables takes the form 
dyj i /dt-Eky z yj i , dy z /dt=E(A + ky 1 , dtp 3 /dl-~w z 


Optimization of the process p° , y° requires the existence of 
such nontrivial constant 0 and vector-function yttj with 


will allow to meet the maximum condition [21_ 

max h(yj(0,y(0,t,p ) - h(y( O >y( t) , t , 
p£ Q 

and transversality condition 




C 2.61 


wvap.Sap.^.p'ctp = l 


C2.7) 


where q = jo, dy z /dt, dy 3 /dt J-| t = t . 

According to C 2.41 the right-hand side of (2.71 is equal 
to zero. 

Hamiltonian operator is expressed as 

h = (v 0 +y 3 )p ~ yE\c\/23Bt z +(A + fcy^y 2 ] + y z y 3 C2.81 


Integration of the system (2.51 combined with conditions 
(2 41, (2.71 enables one to define behavior of Hamiltonian h m 
terms ' of linear function of p with coefficient $ = V 0 + V 3 - 


3. Cylindrical indenter with a conic head of the height 6. 
Linear relation H = <pCp 0. 

The system of differential equations describing the process 
is divided into two parts: 

r -E [U/OB t z + (4 + fey 1 !}/ 2 ], t < 6 

dy /dt - | _ E [U/2 - )B , 5 2 + (A + My‘ly 2 ], t > 6 , (3 1} 


dy z /dt = •{ 


f y 3 , t < 6 


L y 3 - 6pC t-63, t > 6 , 


dy 3 /dt 


■{ 


pCO, t < 6 

pCt D - pCt-6), t > 6 , 


h 


Hamiltonian operator takes the form 
= C^+^pCO " +C4 + My l V] 

h = Cw + W )p(t) - C6ip ~ y^)pet-6) 

T o r 3 J 2 3 

- y^E](i/2)B6 z +C4 + My l 3y 2 ] > 


+ y 3 , for t < 6 ; 

+ v 2 y 3 " 

(3.21 


for t > 6 


The equations for the conjugate variables is expressed as 
.(2 51 y t£[u,t ]. It is assumed that the value y o is such that 

L 
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pondition L k > 6 is satisfied automatically. 

4. Cylindrical indenter with a flat end- face. 

Nonlinear relation H - <pCp). 

The equation of motion for indenter and Hamiltonian are writ- 
ten, respectively as 


dy l /dt = -E i<pdp) + py 1 ] , C4.1) 

h = W Q P ~ W^hpCp) + py 1 ] (4.3) 

Using differential equation for conjugate variables and 
transversal ity condition gives: 

h = C4. 3) 

$ x - i + y 1 fZ t = fZ y Z = exp C-E f k pdrJ C4.4) 

t 


For a piecewise-lmear relation <pC p) (see fig. 2) Hamiltonian 
is transformed to a piecewise-1 inear function p with the slope $ 


-1 + 

$ = < 

~t + 

For nonlinear relation 


(Zp k /(B+Ap k JXA+y l pQlp^pJ 

C Zp k /C B i +Dp k )X D+ y 1 ) , pSTp # , p z J 

<pC pD=B+Ap n CA>0,n>0) Hamiltonian 


(4.5) 
reduces to 


h = pSOpO + $ = t + fZQ y l + Ap n / n = p k /(B+Ap k D (4.6) 

The conditions assumed for existence of continuous solutions 
may be expressed as 

= § + A $ Cp°X~ l = 0 

rs i rN n * 


C dh/dp') 

(d z h/dp z ‘) 

Following (4.7) one gets: 


p=p° 1 n 2 


= A nCn-lX Cp°X 2 < 0 

p=p° n 2 


Differentiat 

dp°/dt = -G(p°X[ 


p° = [-§ /(A $ )] m , 

Jr 1 n 2 J 
ion of (4.9), using 


m. = 1/Cn-l) 
using (4.1) gives 

t + B/(ACl-nX Cp°X n J, 


(4.7) 

(4.8) 


(4.9) 

G = E/n (4. 10) 


RESULTS 

Without going into details we shall examine some qualitative 
results obtained for a number of special cases. 

1. Linear relation H = pCp 0. 

From the analysis of the system (2. 1) , (2. 5) , (2. 8) we can draw 
the following qualitative conclusions : 

(i) There is an interval Q ,i ] in which the optimal functi-J 
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on has the form of p° = p , i.e. the rear layer should be made of 

more light and less hardness materials; 

CiD The function p° = p is optimal within the entire 

interval tt Q ,t k h if Bt 0; 

Ciii) In the case of B < 0 the structure of target should be 
double layer with the front lay being made of hard and heavy mate- 
rial. In this case a relay-type control is realized. 

Similar results have been obtained for the case (3.13,C3.23. 
A general phase diagram of the optimal structure is shown in 

fl9 ' Here v is the velocity at which L k = 6, where 6 the height 

of the indenter head. It is seen that the structure largely de- 
pends on the parameter B. If the specific hardness Q = H/p is 
assumed the measure of material quality, then, according to (2.21, 
the condition B<0 identifies the maximum quality of heavy material 
<p ,H >. In this case a double layer plate is an optimal structure 

for a target. Contrary, when v > v q an optimal structure may be 

represented by a homogeneous plate made from a light material. 


2. Nonlinear relation H = <pCp 0. 


CiD For a piece-wise linear function HCp) as plotted in 
fig. 2 the problem is solved for three different materials. The 
"phase-diagram" of the optimal structure is shown in fig. 4, where 
a = -Cp /p 3Cp -p # 3/Cp # -p i 3, X = (L k /(53 2 -i. As it follows from the 

observable scale effect, the optimal structure for a given set of 
materials CB,B 3 depends on the relation L^/6. 

C i i 3 For nonlinear dependence HCp3 the function of the opti- 
mal control may not include discontinuities. Inequality C4.83 is 
valid for n <1 in the neighborhood of 0<t # <t k . The 

procedure of qualitative estimating the type of solution to the 
equation (4.103 may be as follows. The first approximation (the 
expression in square brackets in (4.103_is constant3 follows from 
dp°/dt = -6Cp°D 2 as Cp°D 1 = p GCt-t } . This solution 


is found to be exact for B =0. The second approximation has a more 
complicated form. The results are shown in fig. 5. The position 
(coordinate3 of the point t 0 is calculated numerically and may 


coincide with the starting point of the process. 

The results of present investigation allow to make a prompt 
qualitative estimation of the optimal target structure. The best 
ratio of layers in a double-layer target may be calculated 
numerically Dy solving the equation of motion for indenter. 
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